Full counting-statistics in a single-electron transistor in the presence of 
nonequilibrium quantum fluctuations of charge 
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Using the Schwinger-Keldysh approach and the drone (Majorana) fermion representation, we 
evaluate current distribution in a single-electron transistor in a regime where a total tunnel resistance 
is small. Nonequilibrium quantum fluctuations of charge induce a lifetime broadening of charge- 
state level, as well as renormalization of system parameters. We find that the lifetime broadening 
effect may suppress the probability of relatively large current. 
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Out of equilibrium, properties of carrier fluctuations 
at zero temperature, T = 0, is beyond the scope of 
the 'fluctuation-dissipation theorem'. They have been 
a highlight of the mesoscopic quantum transport. Re- 
cently, the 'full counting statistics' (FCS) has pervaded 
as a powerful concept for a description of nonequilibrium 
current fluctuations 0. The FCS is attractive because 
it would promote our understanding on nonequilibrium 
strongly-correlated mesoscopic systems. However only 
several attempts have been done for now [3, La- Ll • 

A basic example of the strongly-correlated mesoscopic 
systems is a 'single-electron transistor' (SET), a metallic 
island coupled to left, right and gate electrodes with small 
total capacitance Cs = Cl+Cr+Cg- The island can also 
exchange electrons with left (drain) and right (source) 
electrodes through tunnel junctions with resistance i?L 
and i?R. However, the charging energy Ec = e 2 /(2Cs) 
of a single electron can exceed source-drain bias volt- 
age V (eV > 0) and insulate the island [Coulomb block- 
ade (CB)]. One can reduce energy for impinging an elec- 
tron into a charge-neutral island Ao — Ec(l — 2Qa/e) by 
gate-induced charge Qg- Then electrons start to tun- 
nel through the island one by one [sequential tunneling 
(ST)], when Ao reaches a condition /ir<Ao</il, where 
ML/R = K L /Rey [kl/r = ±C , r/l(C l + Cr)~ 1 ] is a voltage 
drop between the L/R electrode and the island. 

The physics characteristic for strong electron- 
correlations emerges when a total resistance i?x = 
RlRr/(Rl + Rr) reduces to the order of the resistance 
quantum R^ = h/e 2 and quantum fluctuations of charge 
are promoted || . The prime result is renormalization of 
the charging energy and the conductance. In the weak 
tunneling regime where the dimensionless conductance 
a = i?K/(47r 2 i?x) is smaller than unity ao<Cl, the scal- 
ing analysis predicted that the renormalization factor zq 
logarithmically depends on cutoff energy A, which is T 
or eV (h=kB = l), as z = l/{l + 2a ln(£ c /A)} @. 

For SET-like geometry, the FCS was discussed for 
an open quantum dot coupled to electrodes by single- 
channel point contacts Q and for a quantum dot in the 
Kondo regime 0. However, for the metallic island, the 
FCS was obtained only in the limit of cvo — * 7] . The 



FCS in the regime where the total resistance is small is 
still unknown. 

In this paper, we evaluate the FCS of SET in the 
weak tunneling regime non-perturbatively using the 
Schwinger-Keldysh approach and the 'drone' (Majorana) 
fermion representation ||. We will show out of equilib- 
rium, strong quantum charge-fluctuations induce a life- 
time broadening of charge-state level and suppress the 
probability for relatively large current to flow during a 
measurement. Hereafter, we assume the energy and spin 
relaxation times in each region are fast enough and elec- 
trons obey the Fermi distribution f + {u) — l/(e W//T + l). 

FCS in the weak tunneling regime. - In this regime, 
inverse of the RC time 1/r — AitciqEc = 1/{RtC^) is 
smaller than Ec- It ensures that charge-state levels are 
well resolved and the low-energy physics is dominated by 
charge states with charge neutrality and with one excess 
electron. Then a Hamiltonian is mapped onto the 'mul- 
tichannel anisotropic Kondo model' [2j , by introducing a 
spin-1/2 operator a acting on the two charge states, 
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where a r kn annihilates an electron with wave vector k in 
the left or right electrode (r=L,R) or the island (r=I). We 
assumed the tunneling matrix element T r is independent 
of k and the index of channels for transverse mode and 
spin n. The microscopic parameters and the resistance 
of the junction r are related with each other as Rk/Ri = 
(27r) 2 Af c h|T r | 2 pi/? r , where A c h is the number of channels. 
The electron density of states in each region p r is assumed 
to be constant. 

In the nonequilibrium state, the logarithmic renormal- 
ization of the system parameters can be treated by the 
'resonant tunneling approximation' 9] but only for an 
average current. A technical difficulty arises because the 
spin-1/2 operator prevents to utilize Wick's theorem. A 
convenient device to avoid the problem might be the Ma- 
jorana fermion representation [lfj, a_| 



; =2c t c-l 
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where <f> = <$ + d is a Majorana fermion and c and d are 
Dirac fermions. The device enables one to utilize the 
fermionic Schwinger-Keldysh approach 0, 0] . 

A central quantity of the technique is the generating 
functional of connected Green function (GF) 



gen< 
EI 



W[(p] = -ihi J D[a* kn ,a lkn ,c* ,c,d* ,d] expfi J 



dtL{t) ) , 



where C is the Lagrangian for H and we introduced six 
Grassmann variables. The closed time-path C advances 
from t = — oo to oo and returns to — oo. Then it connects 
to the imaginary time path and closes at t= — oo — i/T. 
We introduce an auxiliary source field ip(t), the phase 
of the tunneling matrix element T r — > T T e iKl ' p ( t \ and 
assume a part for tunneling is switched on adiabati- 
cally. Now we can define fields on forward and backward 
branches ip+(t) and tp_(t), but only the center-of-mass 
coordinate <p c (t) = {(p+(t) + ip- (t)}/2 has physical mean- 
ing and is fixed as eVt. 

In the scheme of FCS, during measurement time to, the 
relative coordinate <fA{t) = — is switched on 

and is fixed as a constant A called the 'counting field' [l2j |: 



W(A) - iW[tp}\ iPcW=eVt ^ ip& ( t ) = x8(t o /2+t)0(t o /2-t), (3) 

where 6{t) is the step function. Then the generating func- 
tional is reduced to the cummulant generating function 
(CGF) for the number of transmitted electrons q during 
a measurement, W(A) = E^o ))( iA )"M The dis- 
tribution of q (equivalently current I=eq/to) is given by 
the inverse Fourier transformation, 
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3 W(A*)-i(t J/e)A* 
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The righthand side is a saddle point solution and A* is 
obtained by solving I — ~iedxW(X*)/t 0. The approx- 
imation is valid for long measurement time since W is 
proportional to to as we will mention later. 

We proceed following Ref. 0], where we developed 
a conserving approximation for second cummulants, i.e. 
'noise', based on the Schwinger-Keldysh approach and 
the Majorana fermion representation. We trace out elec- 
tron degrees of freedom and obtain an effective action for 
c and d fields S x = S c h+S x consisting of parts describing 
the charging energy S c h and the tunneling S x 0] as 

S ch = f dt{c(t)*(id t -A )c(t)+id(t)*d t d(t)}, (5) 
Jc 

= - fdtdt' c*(t)<p(t) a x (t,t') <j)(t')c(t') + 0(T r 4 ). (6) 
Jc 

The term 0(T^) describes a phase coherent propagation 
of an electron between the left and right electrodes and 
can be neglected for N c h ^> 1. Here a x = a^ + a^ is a 



particle-hole (p-h) GF describing tunneling of an electron 
from one of the electrodes into the island. In the Keldysh 
space, the p-h GF is expressed as a 2 x 2 matrix rotated 
by A r — K r X along the a;- axis as 



a x (u;) = U{ a r M U K , a T (w) = (j 



af(uj) 



(7) 



(2) where U\ t = exp(— iA r Ti/2) and [Ti]y = 1 



■Sij. Re- 
tarded and advanced components are given by af~(oj) = 
af(u))* = -i7rao {(jJ—Ht) E^/{(oj-fj, T ) 2 +E^}, where a = 
i?K/(47r 2 i? r ) and we introduced the Lorentzian cutoff to 
regularize the ultraviolet divergence. A Keldysh compo- 
nent (uj) = 2a.^(uj){n + (ui — n~(uj — contains 
the Bose distribution function n + (uj) = l/(e"/ T — 1) and 
nT = l+n + . 

We observe that only the Dirac fermion c may carry 
physical meaning and its free retarded GF g^(uj) = l/(w+ 
irj — Ao) describes dynamics of charge-state excitation {jj 
is a positive infinitesimal). Then the self-energy Yj X = 
£l + £r would account for quantum charge-fluctuations 
caused by tunneling. The first order expansion in ao is 
Y, x (t,t') = -ig<f,(t',t)a x (t,f), where g4> is the d-field GF. 
Keldysh and retarded components of the self-energy for 
zero counting-filed S r (w) are given by 



EfH = 2afH, E*(w) = du' 



(2tt)(u> +irj -w')' 



(8) 



The integration results in a logarithm as S^(w) ea 
a ln{2E c /(eV)}w - iY/2 for \w\ < eV (T = and a 
symmetric case, i.e. _Rl = -Rr and Cl=Cr). The imag- 
inary part describes the lifetime broadening of charge- 
state level caused by dissipative current: T is the sum 
of tunneling rates r = riL + rLi + rrR + rRj, where r rI / Ir = 
(A — /x r )n ± (A — fi r ) / (e 2 R r ) is the tunneling rate into (out 
of) the island through the junction r calculated within 
Fermi's golden rule. 

We can perform a systematic perturbative expansion 
in ao for the CGF using diagrams 0|. For exam- 
ple, the first order expansion is given by WW(A) = 
—J c dtdt'g c (t,t')'E x (t',t). We proceed by projecting the 
time on C onto the real axis. For enough long mea- 
surement time, we can approximate an integral St {u>) = 



t /2 

I-t /2 e~ luJt /(2ir) and its square normalized by to, 
2ir5t (u) 2 /to, as the J-function S(lu). Especially, the lat- 
ter equation ensures that W is proportional to- By per- 
forming the Fourier transformation, we obtain 



W [1] (X)^-t J <iwTr{s c (w)Ti£*(u;)"ri}/(27r) 

=*° E r=L . R { p - r ^ elAr -i)+^+rir(c- lA - -1)}, (9) 

where we utilized the expression of Keldysh component 
for free c-field GF g^(cu) = 2ilmgf (w)(P_ - P+). It 
contains initial distribution probabilities of charge state 
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P± = l/(e ±A °/ T + l). Equation © is the sum of CGFs 
of Poissonian distribution describing tunneling of an elec- 
tron into and out of the island. Here, we observe that the 
naive first order expansion Eq. © is insufficient: First, 
it includes the initial distribution probabilities, which 
should be replaced with stationary distribution probabil- 
ities. Second, the CGF should depend only on Al— A r = A 
from the charge conservation . Moreover, the first or- 
der expansion is obviously insufficient to reproduce the 
logarithmic behavior of the scaling analysis. 

Above flaws can be removed by accounting for a sub- 
class of diagrams. Practically, we simply sum up the 
geometric series in (g c TxTi^Ti), which contains leading 
logarithms, i.e. powers of ao ln{2£'c/(eV)}: 



2tt 



W(A) ^^-fdcuTr ln[$ c (w) -n^^n 



dLu\n{l + T"(Lu)f + (Lu-^)f-(Lu- m )(c tX -l) 



+T^a ) )/V- m )/>- / i L )(e- jA -l)} + Ofa),(10) 

where /~ = 1 — / + and we removed a constant to ful- 
fill the normalization condition W(0) = 0. Now P± is 
absorbed in 0(r]) and only A appears in the expression. 
Equation (|10|) looks similar to the Levitov-Lesovik for- 
mula |16(. But for our case the effective transmission 
probability T f {uj) = -a£(u)ag{u)/\u - A - £f H| 2 
accounts for strong quantum charge-fluctuations. 

In fact, Eq. (|10|l is reduction of an approximate 
Keldysh generating functional Eq. (25) in Ref. to 
the CGF. Therefore, first and second cummulants, (I) = 
e((Sq))/to and Su = 2e 2 ((Sq 2 )) /t , reproduce an average 
current and a zero- frequency noise in Ref. |13| . Espe- 
cially, for the average current, our approximation and the 
'resonant tunneling approximation' [2j are equivalent. 

Limiting cases. - In the limit of ao — > 0, we confirm 
that Eq. i|10|) reproduces the 'orthodox' theory 0: 

W«(AM r(\/£(A)-l)/2, (11) 
D(A) = 1 + 4 M^(e lA -l) + 4 %^(e- iA -l). (12) 

It is noticed that though is proportional to ao, it is 
different from the naive first order expansion WW. The 
second order expansion in ao reads, 

W (2) (A) = W cot (A) + a Ao {RcSf(Ao)W (1) (A)}. (13) 

The second term provides the renormalization of the sys- 
tem parameters up to first order in ao |17|. The first 
term is the CGF of a bidirectional Poissonian process 



W cot (A) = £ {7 + (e 4A - 1) + 7" (c~ - 1)} 



(14) 



governed by the cotunneling rate, 7 ± = 2n a}j fdcu(ui— 
/iL)(w-^ R )n ± (w-/i L )n T (w-^ R )Re{(w+^-A ) _2 }. The 
first term is relevant in the CB regime and is consistent 



with the FCS theory of quasiparticle tunneling 18]. Since 
a tunneling quasiparticle is an electron tunneling from 
the left electrode to the right electrode with a p-h pair 
left in the island, the tunneling rate is not proportional 
to |T L/R | 2 but |T L T R | 2 . 

We will discuss the saddle point solution of the cur- 
rent distribution P, Eq. Solid lines in Fig. QJa) 
show the logarithm of nonequilibrium distribution for 
the symmetric case at T — 0. The conductance is 
small and thus and >V cot well approximate CGFs 
deep in the ST regime and in the CB regime, respec- 
tively. At A = 0, the CGF is WW(A) » 2q(e lX / 2 -l), 
where eq/t Q = G Q V/2 and 1/G = Rl + Rr- The fac- 
tor e lA / 2 results in a sub-Poissonian value of the Fano 
factor Sjj/(2e(I)) « 1/2 indicating that tunneling pro- 
cesses are correlated. The meaning of the 'correlated 
process' can be understood from an explicit form of dis- 
tribution P(q) =E^, gR =o p p(lh)Pp(qR) ^( 9L + 9R )/2 ob- 
tained by simply performing the inverse Fourier trans- 
formation Eq. (0J without the saddle point approxima- 
tion. For symmetric case, numbers of transmitted elec- 
trons through the junctions L and R, and g R , follow 
the same Poissonian distribution Pp(q) —q q e~^/ql . The 
Kronecker delta implies that and <7 R are not indepen- 
dent variables. 

As Ao increases and approaches a threshold value Ao = 
eV72, the left junction becomes dominant W^^A) ss 
^oIli (e* A — 1) and the Poissonian value is approached 
Sn I \2e(I)) « 1. The Fano factor continues to be unity 
in the CB regime, because the CGF is also that of the 
Poissonian distribution W cot . However, since tunnel- 
ing quasiparticlcs are different in two regimes, near the 
threshold, there is a regime where neither the 'orthodox' 
theory W^ 1 ' nor the 'cotunneling theory' W cot does not 
work. The plot for Ao/eT^ = 0.475 in Fig.^a) is such an 
example. The 'orthodox' theory [dashed line] cannot fit 
the solid line. 

Renormalization and lifetime broadening effects. 
As the conductance increases, the quantum charge- 
fluctuations are promoted. However, when the life- 
time broadening is negligible, z T «C eV [A = 
max(|z A |,27rT, |eV|/2)], the 'orthodox' theory 
with renormalized parameters z^ao and zqAo is valid. 
Figure QIb-1), the current distribution at Ao = and a 
small bias voltage, indicates the renormalization effect: 
Since zq decreases with increasing ao, the mean value, 
i.e. a peak position, shifts leftwards. The above scenario 
could fail in a regime T K = E c e- 1/{2ao) /(2tt) > A where 
leading logarithms would be not enough. 

If we replot Fig. ^b-1) with vertical and horizontal 
axes normalized not by GqV/2 but by (I), three curves 
almost overlap (not shown), i.e. the renormalization fac- 
tor is normalized away. A qualitative change occurs when 
the bias voltage increases and the lifetime broadening ef- 
fect is enhanced [Fig. njb-2)]. We observe a shrinking 
of the distribution with increasing ao- The shrinking 
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FIG. 1: The zero-temperature current distribution for J?l = 
Rb. and Cl = Cr [(a), (b-1,2)]. (a) Solid lines are plots for 
various excitation energies for small conductance. A dashed 
line is a result of the 'orthodox' theory at Ao/(eV) =0.475. 
Plots at Ao = for various conductances with small (b-1) and 
large (b-2) bias voltages. Axes are normalized by the average 
current {/) (b-2) and that of the 'orthodox' theory GoV/2 (b- 
1). (c) The skewness at eV = 0.1Ec for various conductances. 

is compatible with the suppression of the Fano factor 
predicted in Ref. 01 . The FCS analysis gives more in- 
formation beyond the noise analysis. From the current 
distribution, we can see the probability of relatively large 
current is suppressed. 

Let us discuss the lifetime broadening effect for sym- 
metric case more quantitatively. For eV 3> 7k, the 
real part of £,? can be neglected and only the imagi- 
nary part describing the lifetime broadening is impor- 
tant, (u) « -ina eV. Then the CGF at A = is 
estimated as 

W(A) » 2q{(e lX ' 2 -l) - 2a (e iA -l) 

+ 7r 2 a2(e m /2_e 4A /2)/ 2 + 0( a 3)} ; (15) 

and the 'Fano factor' reads {(8q n )) / {(8q)) = 2 1 -"{1 - 
4ao(2™~ 1 — 1) + O(o!o)}- We can see higher cummulants 
are suppressed as ao increases. Figure Qfc) shows the 
skewness ((Sq 3 )) as a function of Ao. A double-peak struc- 
ture growing with increasing the conductance is obtained. 
Such a characteristic structure might be measurable for 
a present-day experiment [l9j |. 

Summing up, we evaluated the FCS of SET. The 
strong nonequilibrium quantum charge-fluctuations are 
accounted for by summing up a subclass of diagrams in- 
cluding leading logarithms. Our approximation repro- 
duces the 'orthodox' theory in the limit of «o — * and 
gives the bidirectional Poissonian distribution governed 
by the cotunneling rate in the CB regime. Nonequilib- 
rium quantum charge-fluctuations induce not only the 



logarithmic renormalization of the conductance and the 
charging energy but also the lifetime broadening of the 
charge-state level. It may suppress the probability of rel- 
atively large current and higher cummulants of current 
fluctuations. 
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